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The half-filled ground state of the Hubbard model on the hypercubic lattice in D dimen- 
sions is studied by the Kondo-lattice theory, which is none other than the 1/D expansion 
theory, but within the constrained Hilbert subspace where no symmetry is allowed to be 
broken. A gap can open in the single-particle excitation spectrum if and only if the residual 
entropy or entropy at T — +0 K is nonzero. The Mott insulator with no symmetry broken, if 
it is possible, is characterized by nonzero residual entropy or nonzero entropy at T = +0 K. 
This conclusion is consistent with Brinkman and Rice's theory and the dynamical mean-field 
theory. According to the well-known argument based on the Bethe-ansatz solution, on the 
other hand, the half-filled ground state in one dimension is the Mott insulator although 
its residual entropy per unit cell is vanishing in the thermodynamic limit. Two possible 
explanations are given for the contradiction between the present paper and the well-known 
argument. 
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§1. Introduction 



Mott localization or the Mott insulator is a longstanding crucial issue. •' Whether 
*^J | magnetism is itinerant-electron or local-moment magnetism is directly related to the 

physics of the Mott insulator. High-temperature (high-T c ) superconductivity occurs 
in the vicinity of the phase of a local-moment type of antiferromagnetic insulator™ 
Thus, it is anticipated that the mechanism of high-T c superconductivity is closely 
related to the physics of the Mott insulator I2l«»*u 

For example, we consider the Hubbard model at the absolute zero Kelvin. Ac- 
cording to Hubbard's theoryp''''^ 1 when the on-site repulsion U is as large as or larger 
than the bandwidth W, the band splits into the upper and lower Hubbard bands, 
(T) ! i.e., the Hubbard gap opens. According to Gutzwiller's theory,^™ '^ when U > W 

^f-) | and electrons are almost half-filled, electrons are renormalized into heavy quasi- 

particles. According to a slave-boson theorypB a midband appears at the chemical 
potential between the upper and lower Hubbard bands; the midband is none other 
than a band of heavy quasi-particles in Gutzwiller's theory, so that it is called the 
Gutzwiller band. According to Brinkman and Rice's theory,^ which is based on 
^ Gutzwiller's theory, if electrons are exactly half-filled, a metal-insulator transition 

occurs at a critical U c ~ W . In an insulating phase at U > U c , the Gutzwiller 
band disappears and the Hubbard gap becomes a complete gap. According to these 
theories, the ground state is a metal for any U unless electrons are exactly half- filled 
and, therefore, the ground state can be the Mott insulator if and only if electrons 
are exactly half- filled and the onsite U is large enough. 



*' E-mail: fohkawa@phys.sci.hokudai.ac.jp 
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According to Brinkman and Rice's theory,^ the specific-heat coefficient diverges 
as the onsite U approaches the critical U c from the metallic or nongapped phase; 
and no symmetry is broken in the insulating or gapped phase. A combination of 
these two properties means or at least implies that the residual entropy is nonzero 
in the gapped phase, although it is unlikely that the third law of thermodynamics 
is broken. Whether the residual entropy is zero or nonzero in the gapped phase is a 
crucial issue to be elucidated. 

Either Hubbard's or Gutzwiller's theory is under single-site approximation (SSA). 
If all the single-site terms are rigorously considered by a theory, the theory is the 
supreme SSA (S 3 A) theory. In the limit of D — > +00, where D is the spatial di- 
mensionality, multisite terms vanish except for certain types of conventional Weiss 
mean fieldri which is a multisite effect. Thus, the S 3 A is rigorous in the limit of 
D — > +00,^'^'^'™ but within the constrained Hilbert subspace where no sym- 
metry is allowed to be broken. A problem of rigorously calculating all the single-site 
terms under S 3 A can be substituted for a self-consistent Kondo problem 113 ED .LOS M3 
in which the Anderson model to be solved has to be self-consistently determined with 
the single-site terms. A dynamical mean- field theory^ •^ (DMFT) or a dynami- 
cal coherent-potential approximation (DCPA) theory^ 1 is none other than the S 3 A 
theory. 

Under DMFT or S 3 A, a metal-insulator transition occurs at U as large as 
Brinkman and Rice's [/ C J22IJ23)J23 Since hysteresis appears, it is a discontinuous 
transition between the metallic and insulating phases at T = OK. Since no jump 
occurs in the internal energy of Brinkman and Rice's theory, a plausible explanation 
for the hysteresis is that the residual entropy jumps at the transition. It is interesting 
to elucidate which is discontinuous at the transition point and is responsible for the 
hysteresis, the residual entropy, the internal energy, or some thing else. 

Even if multisite terms are approximately or rigorously considered beyond S 3 A, 
a problem of rigorously calculating all the single-site terms can also be substituted 
for another self-consistent Kondo problem™ *B& ,LL3),|26) - m ^-^j^ th e Anderson model 
to be solved has to be self-consistently determined with multisite terms as well as 
the single-site terms. Since the Anderson model is an effective Hamiltonian for the 
Kondo problem, a theory for treating either the self-consistent Kondo problem under 
or beyond S 3 A is none other than a Kondo-lattice theory IIB>II2^LL9{,126J> jf multisite 
terms are considered beyond S 3 A, the Kondo-lattice theory is none other than a 
1/.D expansion theory. In the present paper, we list all the possible types of gapped 
phase of the Hubbard model under and beyond S 3 A and then we examine on the 
basis of the Kondo-lattice or 1/D expansion theory whether the residual entropy 
is zero or nonzero in each of them; we consider only possible phases within the 
constrained Hilbert subspace in which no symmetry is allowed to be broken, so that 
the possibility of the Mott insulator with symmetry broken or one with a hidden 
order is out of scope. 



*' In a single-band model such as the Hubbard model, three types of conventional Weiss mean 
field are of the zeroth order in 1/D: spin density wave (SDW) or magnetism, charge density wave 
(CDW), and Bardeen-Cooper-Schrieffer (BCS) or isotopic s-wave superconductivity. In a multiband 
model, an orbital order in the charge or spin channel is also of the zeroth order in 1/D. 
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The present paper is organized as follows: Preliminaries are given in £}2j In £}3l 
nongapped and gapped phases under S 3 A are studied by the Kondo-lattice theory in 
order to confirm that, under S 3 A, the Kondo-lattice theory is equivalent to DMFT. 
In §U only possible complete-gapped phases beyond S 3 A are studied by the Kondo- 
lattice theory. Discussion is given in $5j Conclusion is given in £j6j Anomalies 
of possible gapped phases are studied in AppendixlAl The entropy at T = +0 K of 
possible complete-gapped phases is studied in AppendixlBl Possible gapped phases in 
a nonlattice or continuous model, such as the Tomonaga-Luttinger mo del l2ZJJ2a)J22) 
are studied in AppendixICl 

§2. Preliminaries 

2.1. Model 

We consider the Hubbard model on a hypercubic lattice in D dimensions: 

U = -^ tijd\ a d j(J + U^2 «it n 4> ( 2-1 ) 

ijcr i 

where n^ = d\ a di a and the other notations are conventional. The periodic boundary 
condition is assumed. When U = 0, the dispersion relation of an electron is given by 



E(k) = -j^t ije ik < Xi ~ x ^ , (2-2) 



L 
v 

where k = (k\, k,2, ■ ■ ■ , kr>) is the wave number of the electron, Xi = (i\, t2, ■ ■ ■ , i_o)a, 
with a being the lattice constant, is the position of the ith unit cell, and L is the 
number of unit cells. The hypercubic lattice is decomposed into two sublattices: one 
composed of unit cells with X^=i *&> being even and the other composed of unit cells 
with J] , i v being odd. For simplicity, we assume that the transfer integral tij is 
nonzero only between different sublattices, i.e., the Hubbard model is on a bipartite 
lattice. The density of states per unit cell is defined by 

Po(e) = ij>[e-tf(fc)], (2-3) 

k 

which is called the bare density of states in this paper. The band center of E(k) or 
po(e) is denoted by e^; e^ = for the model of this paper. A hypersurface in the 
wave- number space, k on which is denoted as kp, is defined by 

E(k F ) = e d = 0. (2-4) 

Since E(k) is symmetrical with respect to the hypersurface, po{e) is symmetrical 
with respect to e^ = 0: po(^) = Po(— £)• We assume that the chemical potential p, 
is exactly given by p = U/2 and that the Hubbard model is exactly half-filled. In 
general, the particle-hole symmetry exists in a half-filled model on a bipartite lattice. 
The particle-hole symmetry is crucial in this study. 



4 F. J. Ohkawa 

The transfer integral tij has to include the dimensional factor D in order for the 
effective bandwidth of E(k) to be finite for any DjMHMSBMb p r example, if tij is 
nonzero only between nearest neighbors, tij = t/yD between them. It follows that 

t D 
E{k) = -1—j= Yl cos (^ a )- ( 2 ' 5 ) 

The absolute bandwidth of E(k) is given by W = 4|t|-/D. Since (1/L) ]T fc E ( k ) = ° 
and {l/L)J2 k E 2 (k) = 2t 2 , the effective bandwidth of E(k) is O (\t\) for any D. If 
t^ is of finite range, either the absolute or effective bandwidth is finite for any finite 
D. We call any model of this type for a finite D a finite bandwidth model. 
If Uj is of infinitely long range, it is possible that E(k) is given by 

£(fc)=E%W*), (2-6) 



, i D 



where 



) tan (x - tt/2) , < x < vr, 
^W = S / , x (2-7) 

) - tan (x + tt/2) , -vr<x<0. 

The bare density of states is of a Lorentzian type for any D: 

J-i " T K-oo 

We call this model an infinite bandwidth model. 

2.2. Kondo-lattice theory 

We constrain the Hilbert space within the subspace where no symmetry is al- 
lowed to be broken. We first assume that L is large but still finite and that the 
temperature is nonzero such that T > OK; the thermodynamic limit of L — > +oo, if 
necessary, followed by the low-temperature limit of T — > K is taken last. 

The site representation is convenient for formulating the Kondo-lattice theory. 
The Green function and self-energy in the grand canonical ensemble are denoted 
by Rij- ff (isi) and Uij- a (iei), respectively, where E\ = (21 + l)i:k-QT, with I being an 
integer, is a fermionic energy. We consider a skeleton Feynman diagram for Eij- a (iei), 
but with each electron line replaced by a renormalized line of Rij- a (i£i). If only lines 
of U and site-diagonal Ru a (iei) appear in it, it is a single-site diagram; if at least one 
line of site-off-diagonal Ri^j. a (iei) appears in it, it is a multisite diagram. All the 
Feynman diagrams can be classified into the single-site and multisite ones without 
any ambiguity. Thus, 

Zij;a( i£ l) = s i;cr(i£l)Sij + AEij :a {ie{) , (2-9) 

where Ei ;(T (iei) is the single-site self-energy, 6ij is the Kronecker delta, and AEij- a {ie{) 
is the multisite self-energy Since the site-diagonal Rn ;a (iei) and the single-site 
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£i;<r(i£i) do not depend on the site index i, they are simply denoted by R a (iei) 
and S a (isi), respectively, in the following part of this paper: 

The single-site E a {ie{) is of the zeroth or leading order in 1/D, while the multisite 
AUij :a (isi) is of higher order in 1/£)L12JA12)A1£JA1§J Any type of conventional Weiss 
mean field is a multisite effect. If no multisite term is considered, the Kondo-lattice 
theory is reduced to S 3 A, which is also DMFT or DCPA. In this paper, the multisite 
AEij., ff (isi) is considered, but no conventional Weiss mean field is considered. 

If the multisite ASij :(T (i£i) is ignored under S 3 A or if it is approximately or 
rigorously considered beyond S 3 A, the single-site U a {iei) is exactly equal to the 
local self-energy E a (iei) of the Anderson model that is determined by solving the self- 
consistent Kondo problem, as discussed below. The mapping of S a {ie{) to IJ a (iei), 
such as 

E a (iefi = S {iei), (2-10) 

is the key of the Kondo-lattice theory. The Anderson model is defined by 

H = £d ^2 hd ° + ^2 Ec{k)c\ a c ka + —= ^2 ykc\ a d a + V£d\c ka \ + Uh d ^h di , 



L ka 

(2-11) 

where h da = d a d a and the notations are conventional except for the use of tildes. 
The local Green function for localized electrons in the Anderson model is given by 

G*(m) = . — : ~ r. : W . t , (2-12) 

iei + fx-e d - 2J a (iei) - r(iei) 

where jl is the chemical potential, E a (iei) is the local self-energy, and 

1 f +ao Ate) 
r(iei) = - de — ^-, (2-13a) 

vr J-oo m ~ e 



with 



A(e) = jYl \^\ 26 i e + ^ " Ec(k)} , (243b) 

L k 



is the hybridization term with conduction electrons; and A(e) = — ImT(£ + iO). 
Although the Anderson model is characterized by jl — e d , U, E c (k), and Vk, it is 
essentially uniquely characterized only by ji — e d , U, and A(e). In the Kondo-lattice 
theory, not only the local U a (iei) of the Anderson model and the multisite ASij :(T (iei) 
of the Hubbard model, but also the jl — e d , U, and A(e) of the Anderson model have 
to be self-consistently calculated or determined to satisfy Eq. (|2T0p : the Anderson 
model determined in this way is called the mapped Anderson model. 

Even if any arbitrary diagram for U a {iei) of the Hubbard model is considered, 
the very same one for E ff (iei) of the Anderson model is obtained. This one-to- 
one correspondence between the diagrams for E a {ie{) and S a {iei) is unique, and 
there is no lonely diagram in either S a {ie{) or E a {iei). In each pair of diagrams, 
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the interaction and electron lines are U and R a (iei), respectively, for U ff (iei) of 
the Hubbard model, while they are U and G a (iei), respectively, for £ a (iei) of the 
Anderson model. Provided that a couple of 

U = U, (2- 14a) 

and 

G a {iei) = R a {iei), (2-14b) 

are satisfied, the contributions of any pair of diagrams are exactly the same as each 
other; thus, the mapping of Eq. (|2-10p is satisfied. The condition of Eq. (|2-14bp is 
equivalent to a couple of 

fi - e d = /J, - e d , (2- 14c) 

where e^ is the band center and ea = 0, and 

A(e) = lm £ a (e + i0) + l/R a (e + i0) , (2-14d) 

where S ff (e + iO) and R a (e + iO) are the analytical continuations of IJ a (i£i) and 
R a (iE[), respectively, from the upper-half complex plane onto the real axis. Since 
Eqs. (|2-14ap and (|2-14cp are trivial, Eq. (I2-I4dll is a practical mapping condition for 
determining the mapped Anderson modelHI 

If a series of Feynman diagrams, single-site or multisite ones, for the Hubbard 
model are divergent in the thermodynamic limit of L — > +oo followed by T — > K, 
the limit has to be taken at this stage. The mapped Anderson model depends on L 
and T. In the following part of this paper, we consider the mapped Anderson model 
in the limit of L — > +oo, if necessary, followed by T — > OK. 

In the wave- number representation, the self-energy is given by 

S a {ie h k) = S^iisi) + AE a {ie h k), (2-15) 

where Eq. (|2-10p is used for the single-site self-energy, and 



L 



The Green function is given by 



GA,e, fe) = | E e -"<^>iW^> = fe|+ „_ E( t ) _ a(fe „ lt) - P-"> 



The site-diagonal Green function is given by 



Ml£l) ~ I ^ iei + ^-E(k)-K(iei,kY (2 ' 18) 



*' Following the previous paperpS* it can be proved on the basis of Eq. (|2-14dl) that A{e) > 
for any possible self-consistent solution, as anticipated. 
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The density of states is simply given by 

p{e) = Imi? CT (e + iO). (2-19a) 

Because of Eq. (|2-14b|) . it is equal to that of the mapped Anderson model either 
under or beyond S 3 A: 

p(e) = lmR a (e + iO) = ImG a (e + iO). (2-19b) 

7T 7T 

This equality is equivalent to the mapping condition and is crucial in this paper. 

When (j, = U/2, not only the Hubbard model but also the mapped Anderson 
model is half-filled and symmetrical: 

H = fL-€ d = U/2 = U/2, (N)/L = (h t + h i ) = l, (2-20) 

where M = ^ icr ni a and (• • • ) stands for the statistical average. Since the particle- 
hole symmetry exists in any half-filled symmetrical model, various equalities are 
inevitably satisfied: In particular, 

Rei7 CT (+i0) -/i = 0, Ref(+i0) = 0, E(k F ) + ReZ^+iO, fc F ) - p = 0, (2-21) 

are crucial, where kp is defined by Eq. (|2-4|) or E(k-p) = 0. 

It is easy to see from Eqs. §FTH) . (pdHll . (pT9l> . and (l2 : 2T]l that, unless Jm[S ff (e+ 
iO)+r(s+iO)\ is divergent at e = 0, then p(0) > and not p(0) = in the Anderson 
model and that, unless ImU a (e + iO, fcp) is divergent at e = 0, then p(0) > and not 
p(0) = in the Hubbard model. A couple of these features are also crucial in this 
paper. 

2.3. Kondo effect 

2.3.1. Normal Fermi liquid in the Anderson model 

The Kondo problem has already been solved EOJ ,[3T|f ,[32j ,[33j ,[34j ,[35j ,[36j j n particular, 
the Bethe-ansatz solution was given to the Anderson model in which A{e) is constant 
as well as to the s-d model in which the density of states of conduction electrons is 
constant W^ <GS)>122I>:SQ} p or exam ple, we consider the Anderson model defined by Eq. 
(|2-lip . The static spin susceptibility of localized electrons is given by 

Xs(T) = (V%Vb lxn( T ) - xn( T )} , ( 2 - 22 ) 

where g and ps are the g factor and Bohr magneton, respectively, of localized elec- 
trons, and 

w(T) = y dr(e^-^)W-^-^rw), (2-23) 

where N = ^2 a fida + J2ka ^ka^ka- Since the most crucial physics of the Kondo effect 
is the quenching of a local moment by conduction electrons J2SDJSDJ22) determining 
whether the Fermi surface of conduction electrons exists is crucial. The Fermi surface 
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is defined by jj, = E c {k). If the Fermi surface exists and | Vfc| 2 is nonzero at least on 
a part of the Fermi surface, it follows that 

< A{e) < +00, (2-24) 

at least for e — 0, so that the ground state is a singlet and a normal Fermi liquid. 
The Kondo temperature 7k is defined bjo 

k B T K = (l/4)g 2 i/ B /xs(0). (2-25) 

If T <C T&, Xs(T) exhibits the Pauli paramagnetism; if T S> Ik, Xs(T) obeys the 
Curie- Weiss law, i.e., localized electrons behave like a local moment or a free spin. 

If ^(0) is neither zero nor infinite and if Ik is nonzero, the ground state is a 
normal Fermi liquid, so that the self-energy of localized electrons is analytical at 
e = and it can be expanded in a way such thaP^'ESU'Eii 

£r(e + t0) =hu+{l-4> x )e- ifoe 2 /A(0) ~ #3 (k B T) 2 /A(Q) + • • • , (2-26) 

for \e\ <C k B T K and T <C T K , where & > 1, 4> 2 > 0, and <fe > 0. When T = OK, 
according to Eqs. ([242]) . (pT9|) . (pF2"T]) . and (f2^6|) . the density of states of localized 
electrons at the chemical potential is simply given by 

p(0) = l/[7rA(0)]. (2-27) 

If A(e) is constant such that A{e) = ^(0), Anderson's compensation theorem is 
validp^ 1 there is no additional contribution of conduction electrons when either V& 
or U is introduced. Then, the expansion coefficients are simply given by^'"^'"^' 

k = xtt(o) M(o)] , h = h = {xn(o) M(o)] } 2 /2- (2-28) 

The specific-heat coefficient of localized electrons is given by 

7 = (2/3)n 2 kUi/ M(0)]. (2-29) 

The Fermi-liquid relation above is consistent with the Bethe-ansatz solution.^* J3SI J39| J40J 

If A(e) depends on e, the compensation theorem is not valid, in general; there 
are additional contributions of conduction electrons. If the compensation theorem 
is valid or not valid, the static spin susceptibility and specific-heat coefficient of 
localized electrons, which include no additional contribution of conduction electrons, 
are given by Eqs. (|2-22p and (|2-29j) . respectively, but with 0i / xtt(0) [^(0)] in 
general. The Wilson ratio is defined by 

w._mw?m, (2.30, 

7 (1/2)^/4 
which is of the order of unity or 0(1), at least, when the onsite U is repulsive. 



*' This definition of Tk is different from Wilson's^ by a numerical factor, even in the s-d limit. 
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If Tk is nonzero, the ground state of the Anderson model is a normal Fermi liquid. 
When T = OK, Imi7 (T (+iO) = and the entropy or residual entropy of localized 
electrons is zero. When OK < T <C Tk, Imi7 o -(+i0) is negatively nonzero and small, 
and the entropy of localized electrons is much smaller than k-Q In 2. When T <C Tk, 
localized electrons never behave like a free spin. When T > Tk, ImZ' o -(+i0) is 
negatively large, the entropy of localized electrons is as large as k-Q In 2, and localized 
electrons behave like a free spin. This crossover as a function of T is crucial.^ 

2.3.2. Possible gapped phases in the Anderson model 

If Eq. (p 7 ^) is not satisfied, it is possible that T K = OK or T K = +0K. If 
Tk = K or Tk = +0 K, the self-energy E ff {e + iO) is singular at e = 0. We examine 
how singular E a {e + iO) is in the case of Tk = K or Tk = +0 K. 

According to Hubbard's theoryp* in the so-called atomic limit of W/U — > 0, the 
self-energy of the Hubbard model is given by 

Since Hubbard's theory is under SSA, it is anticipated on the basis of S 3 A or the 
Kondo-lattice theory that U a (e + iO) given by Eq. ()2-3ip can be approximately used 
for E a (e + iO) of the Anderson model. The atomic limit in the Hubbard model 
corresponds to the limit of Tk — > K or the case of Tk = K in the Anderson 
model. It is anticipated that, if Tk = +0K or Tk = OK, the self-energy E ff (e + iO) 
of the Anderson model has a pole at £ = 0. 

Next, we examine how the self-energy E a (e + i0) of the Anderson model becomes 
singular as Tk ->0K. For example, we consider a function defined by 



1 U 2 f +co dn r/ 2 
E a (e + i0; l0 ) = -U + — -L / d9 C(lho,0) S(e;j,9), (2-32) 

2 4 J +0 70 J +0 



i rr , u 2 r +QO d^_ ri 2 

l+o 7o J+o 
where 70 — &bTk is a parameter, £(x,#) is a functional parameter, and 

-i(7T-0) 



S(e; 7 ,9) ' 



2 cos 9 



-e 



+ 



e — 7e _l ( 7r e ) e — 7e 



-i6 



-e(7 - £ ) - i2£ 2 jsin9 
e 4 - 2e 2 7 2 cos 29 + 7 4 

(2-33) 



Here, we assume that C( x ,6) satisfies C( x ,6) > 0) C( x ~^ +00, 9) = 0, and 

/+00 /"f/2 r+00 /*""/2 

dx d9((x,9)= dx I d9((x,9)x = l. (2-34) 

,-0 J+o J+o J+o 

It is easy to see that S(e; 7, 9) as a function of e is analytical in the upper-half 
complex plane at least for 7 > and < 9 < it/2. It should be noted that 

1 f + °° , 27e 2 sin6' 

de- A ihi ~„ J = 1- 2-35 

£ A - 2e 2 7 2 cos 29 + j A y ' 
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Thus, — (l/7r)ImS'(e; 7 — > 0, 9) is the delta function 5(e) or, at least, it can be treated 
as the delta function, although ImS'(e = 0; 7, 9) = for any nonzero 7. It follows 
that 



1 TT U 2 (_ .. e 2 

e/70^0 2 47^ V 7o 



lim __ E a {e + i0; 70) = -U - -^ ( s\£ + is-i— H ) , (2-36a) 



lim £ CT (e + iO; 7o ) = \u + — ( - - is 3 -% + ■ ■ • ) , (2-36b) 

e/70— >-±oo z 4 ye e 



where 



r+00 rn/2 1 /.+00 /-Tr/2 2 Sm 6* 

Ii= / dx / dffC(ic,^)-Q, S2=/ tic/ d9((x,9) — 3—, (2-37a) 

s 3 = dx d9((x,9)2xsm6. (2-37b) 

i+o i+o 

If 70 ~ /cbTk and Q{x,9) are properly chosen, i7 ff (e + iO;7o) behaves in almost the 
same way as E a (e + i0) given by Eq. (|2-26p for |e| <^ 70 or |e| <C IsbTk- On the other 
hand, it follows from Eqs. f|2-35|l and ()2-36bp that 

lim lim E^e + .0; 70) = \u + ^--^. (2-38) 

Thus, i7 CT (e + iO;7o —5- 0) has a pole or a variation of the pole at e = 0; the pole is 
not simple if 70 is nonzero. It is anticipated that the analytical self-energy becomes 
singular in this way or a similar way as 70 — > or Tk — > K. 

It is anticipated on the basis of the arguments above that there is no essential 
difference between the self-energy characterized by Tk = K and that by Tk = +0 K. 
If the self-energy E a {£ + iO) has a pole or the variation of the pole discussed above 
at e = 0, if the residue of the pole is large enough, if ImZ' ff (e + iO) = for |e| < e g /2 
but £ = 0, and if A(s) = for |e| < e g /2, a gap as large as e g opens in the Anderson 
model0 

To be precise, it has never been fully clarified what singularity the self-energy 
E a (e + iO) of the Anderson model has to have at e = in the case of Tk = K and 
Tk = +0K, i.e., a simple pole, a variation of the pole, or some thing else. However, 
at least, it is certain that, unless ^V(O) is exactly zero or infinite and unless Tk is the 
absolute zero or infinitesimal, the self-energy E a (e + iO) has no type of singularity 
at e = and it can be expanded as in Eq. ()2-26|) , 



Unless A(0) is infinite !**^ a gap can open only if at least ImE ff (£+i0) is divergent 
at s = 0. The divergence of Imi7 CT (e + iO) at e = is possible only if Tk = OK or 
Tk = +0K. If Tk = OK, the residual entropy is nonzero; localized electrons behave 
like a free spin even at T = K. If T = +0 K, the entropy is zero at T = K but is 



*' According to Hubbard's theoryp'" when U/W is so large that a Hubbard gap as large as e g 
opens, E a (e + iO) has a pole at e = 0, the residue of the pole is large enough, and \mE a (e + iO) = 
for |e| < e g /2 but e — 0, as in the Anderson model. 
**' If A(0) is divergent, a gap inevitably opens. 
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nonzero at T = +0K; localized electrons behave like a free spin even at T = +0K. 
The third law of thermodynamics is broken in any gapped phase of the Anderson 
model characterized by 7k = K or Tk = +0 Ko 

2.3.3. Kondo temperature of the mapped Anderson model 

According to the Kondo-lattice theory, the Kondo temperature Tk or &bTk of 
the mapped Anderson model is al so t he energy scale of single-site quantum spin 



fluctuations in the Hubbard model l**1 Thus, we also call the energy scale in the 
Hubbard model the Kondo temperature; it is also denoted by Tk or /cbTk- If the 
Kondo temperature Tk of the mapped Anderson model is the absolute zero or in- 
finitesimal or if the self-energy U a (e+i0) of the mapped Anderson model is divergent 
at e = 0, it can be concluded that at least the single-site part of the residual entropy 
or entropy at T = +0 K is nonzero in the Hubbard model. 

§3. Nongapped and gapped phases under S 3 A 

3.1. Infinite bandwidth model 

We consider under S 3 A the infinite bandwidth model in D > 1 dimensions, in 
which the bare density of states po{E) is given by Eq. (|2-8p . Since po(E) is of a 
Lorentzian type, the problem is exactly the same as that for the Bethe-lattice model 
in D — )■ +oo dimensions™ which is exactly solvable. 

The multisite AE a (iei,k) vanishes under S 3 A. It follows from Eq. (|2T8p that 

/+oo -1 
dEp (E) — -=- — (3- la) 

-oo e + n- E - £ a (e + iO) 

= = — . (3-lb) 

e + fi - E a (e + iO) + i\t oo\ 

It follows from Eq. (|2-14d|) that A[e) = |too|- The self-consistent Kondo problem is 
reduced to a non- self- consistent Kondo problem, to which the Bethe-ansatz solution 
was given W* J23M1IJI2} The ground state is a normal Fermi liquid; Tk is nonzero. The 
density of states is simply given by p(e) = — (l/ir)ImR a (e + iO), where R (7 (e + iO) is 
given by Eq. (|3-lbjl . Since the single-site S a (e + iO) can be expanded as shown in Eq. 
(|2-26p . the density of states in the vicinity of the chemical potential is approximately 
given by that of a Lorentzian type, whose bandwidth is 0(/cbTk). In particular, 
p(0) = po(0) = l/(7r|too|) at T = OK. When T = OK, p(0) is never renormalized by 
U. 



*' One of the expressions of the third law of thermodynamics is that the residual entropy has 
to be zero. Another expression is that the absolute zero Kelvin or T = K can never be reached. 
Since T — K can never be reached, it is possible to extend the first expression to the third one: 
Entropy has to be continuously decreasing and vanishing as T — > K. ff the possibility of an ordered 
ground state in an actual system in three dimensions is considered, the third expression of the third 
law can be extended to the fourth one: Entropy has to be continuously decreasing and vanishing 
faster than linearly inTasT->0K. 

**' Since the mapped Anderson model is determined for a given T, it depends on the given T, in 
general. Thus, the Kondo temperature Tk in the Hubbard model depends on T, in general. 
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3.2. Finite bandwidth model 



When the bare bandwidth W is finite and the onsite U is large enough such that 
U > W, a complete gap opens at T = OK under DMFTpS'GU'IZO 1 which is simply 
S 3 A. In this subsection, we list all the possible types of solution under S 3 A and we 
examine the nature of each of them. We do not specify a precise model for po(E); 
we only assume that p (E) > for \E\ < W/2 and p (E) = for \E\ > W/2. 

The density of states is given by Eq. (|3Tap for the Hubbard model and 

p( £ ) = Im s = , (3-2) 

vr e + /i _i; ff ( e + io)-r(e + iO) 

for the mapped Anderson model. Under S 3 A, E a (s + iO) and r(e + iO) are self- 
consistently calculated or determined in a way such that Eqs. (|3-lap and f|3-2|) are 
equal to each other. According to Eq. (|2-13p . if one of the real and imaginary parts 
of r(e + iO) is finite and continuous at e = £o, the other is also finite and continuous 
at e = eoj an d if one of them is discontinuous or divergent at e = £o, the other is 
divergent or discontinuous at e = eq. In general, 

~ . 1 f +oc , TmEJx + iQ) , x 

S a (E + iO)-p = / dx CT : n ' , (3-3) 

IT J-oo E + lO-X 

for any symmetrical model. According to Eq. (j33fl . the real and imaginary parts of 
Ea(£ + iO) satisfy the very same relation as that for r{e + iO) discussed above. 

First, we consider a possible solution in which HeU a (£ + iO) is finite and con- 
tinuous at £ = 0; ImZ' cr (e + iO) is also finite and continuous at e = 0. Since 
ReZ' cr (+iO) — p = and Re.T(+iO) = because of the particle-hole symmetry, 
it follows from either Eq. (pFTal) or fl£S|) that p(0) > 0. Since E a (e + iO) is fi- 
nite and continuous at £ = 0, according to the mapping condition of Eq. (j214d|) . 
A(e) = (— l/7r)Im-T(e + iO) is nonzero, finite, and continuous at e = 0; Re-T(e + iO) 
is also finite and continuous at e = 0. Then, p{e) > for e ~ 0. Neither possibility, 
Imi7 cr (+i0) < nor Im.E a {+iQ) = 0, can be excluded. According to the correspon- 
dence to the mapped Anderson model, any solution characterized by ImI7 o -(+i0) = 
is one for a normal Fermi liquid at T = OK and any solution characterized by 
ImIJff(-j-iO) < is one for a normal Fermi liquid at T > OK. If ImZf ff (+i0) = 0, it 
follows from Eq. (|3-lap that p(0) = po(0) for the Hubbard model and it follows from 
Eq. (HEp that p(0) = l/[vrZ\(0)] for the Anderson model, as shown by Eq. RFF!\\ ; 
thus, ^(0) = l/[7rpo(0)]. Neither p(0) nor A(0) is renormalized by U. As discussed 
in the following paragraphs, a gap inevitably opens if Rei7 (T (e + iO) is discontinuous 
or divergent at e = 0. Thus, the possible solution examined in this paragraph is 
the only possible one for a nongapped phase under S 3 A, so that it is none other 
than the metallic solution obtained under DMFT for a sufficiently small U, i.e., 
U < ^7 E2J) ,[24J ,[25J) Thjg characterization is never inconsistent with the feature of the 
metallic phase under DMFTJ^'''^'^'' or rather consistent with it; for example, the 
constancy of p(0) = po(0) and A(0) = l/[7rpo(0)] at T = OK as a function of U is 
realized in DMFTpJ^-ES 
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Second, we consider one in which Rei7 -(e + iO) remains finite as e — > ±0 but it 
is discontinuous at e = 0. It follows from Eq. (|3-3p that 

lim ImZUe + iO) = -oo. (3-4) 

Since ImX! cr (e + iO) is nonzero at least for a sufficiently small \e\, if this type of 
solution is possible, no complete gap but only a pseudo-gap or a zerogap opens. 

Third, we consider one in which ReU a (e+iO) as well as Im£ a (e-\-i0) continuously 
diverges as e — > 0: 

lim ReS a (e + iO) = — oo, lim ReS a (e + iO) = +oo, (3 - 5) 

as well as Eq. ()3-4p . If this type of solution is possible, a zerogap opens in it. 

Fourth, we consider one in which KeU a (e + iO) discontinuously diverges at e = 0, 
i.e., ReZ'o- (e + iO) includes the delta function (5(e). In this case, S a (e + i0) cannot be 
analytical in the upper-half plane. No self-consistent solution of this type is possible. 

Last, we consider a possible solution in which ReU a (e + iO) continuously diverges 
ase-> ±0 and \ml! (7 {e + iO) discontinuously diverges at e = 0, i.e., S a (e + iO) has 
a pole at e = such that 

E <r (e + iQ)-(i=-^- + ---. (3-6) 

Then, p(0) = in this type of solution. A necessary condition for a complete gap 
as large as e g to open in the Hubbard model is that S a (e + iO) has a pole at e = 0, 
the residue of the pole is large enough, as discussed below, and ImZ , (T (g + iO) = for 
|e| < e g /2 but e = 0. There is a zero point ze such that it satisfies 

£ + V - z E - Re£ a (e + iO) = 0, (3-7) 

for any finite e but e = 0. Since po(E) > for \E\ < W/2 and po(E) = for 
\E\ > W/2, \e + p- ReIJ a (e + i0)\ > W/2 for any e of |e| < e g /2 has to be satisfied 
for a complete gap as large as e g to open. Thus, the bandwidth W has to be finite 
and the residue Xs has to be as large ao 

\z>\{e 2 g + e g W). (3-8) 

If the bandwidth 2W is finite, a possibility that a complete gap opens can never 
be excluded. This type of solution is the only possible one for a complete-gapped 
phase under S 3 A, so that it is none other than the insulating solution obtained under 
DMFT for a sufficiently large U, i.e., U > py)22IJ2^)J2S) This characterization is not 
inconsistent with the feature of the gapped phase under DMFTp^'^ 1 '^ or rather 
consistent with it, as discussed below. 

In this subsection, we consider an unspecified finite bandwidth model. Any pos- 
sibility of nongapped, zerogapped, and complete-gapped phases cannot be excluded 



*' If the bare bandwidth is infinite such as W = +oo, it is easy to see from Eqs. (|3-lap and (|3-7|) 
that no complete gap but only a zerogap can open. 
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within the qualitative analysis of this section. However, it is certain that ImZ , CT (+iO) 
is divergent in any possible type of gapped phase. Thus, it can be concluded on the 
basis of the analysis in £ )2.3.2l that Tk = K or Tk = +0 K in any gapped phase and 
that the residual entropy or entropy at T = +0 K per unit cell is therefore nonzero 
in it. As regards a complete-gapped phase, in particular, the entropy at T = +0K 
is &b In 2 per unit cell, as studied in AppendixlBl 

In any possible solution of a zerogap or complete-gap phase, inevitably A(0) = 
under S 3 A, as studied in AppendixlAl as well as p(0) = 0. The feature that p(0) = 
and A(Q) = in the possible gapped phase is different from the feature that either 
p(0) = po(0) or A(Q) = 1/[ttpq(0)] is nonzero and constant in the possible nongapped 
phase; therefore, irA(0)p(0) = 1 in it. If a metal-insulator transition occurs at 
T = K or T = +0 K under S 3 A, not only the residual entropy or entropy at 
T = +0K, but also p(0) and A(0) jump at the transition. The jumps of p(0) and 
Z\(0) are consistent with those in the metal-insulator transition in DMFT.'^ 3 -''^''''^' 

§4. Gapped phase beyond S 3 A 

In this section, we consider only possible complete-gapped phases at T = K or 
T = +0 K beyond S 3 A in a finite or infinite bandwidth model. A precise model for 
Po(s) is neither specified. 

The single-site Green function of the Hubbard model is given by R a {e + iO) 
defined by Eq. (|2-18p and the local Green function of the mapped Anderson model 
is given by Eq. (J2-12J) . i.e., by 

G a {e + »0) = -~ , (4-1) 

V ' e + i0-S(e + i0) V ' 

where 

§(e + »0) = [E a {e + *0) - p\ + f(e + tO). (4-2) 

If the multisite AU a (e + i0,k) is normal or anomalous, S ff (e + iQ) and r(e + i0), i.e., 
S(e + i0) of the mapped Anderson model is determined in a way such that R a (e + i0) 
and G a {e + iO) are equal to each other: R a (e + iO) = G a (e + iO), as shown in Eq. 
(|2-14bj) . In all possible solutions, it follows that 

1 f +co o(x) 

R a {e + iO) = — / dx P{ J , (4-3) 

vr J_ 00 e + iO-x 

where p(e) = (—l/-K)lmR CT (e + iO) is the density of states, and 

S{e + i0) = ~L dX e + \o-x - (44) 

We assume that a complete gap as large as e g opens in the p(e) of a self-consistent 
solution: p(e) = for |e| < e g /2. Then, Eq. (|4-3p can be described as 



R (T ( £ + i0) = / dx+ / dx — ^ — . (4-5) 
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Then, for \e\ < e g /2, R a (e+iO) is real and (d/de)R a (e+iO) < 0. Since ReR a (+iO) = 
and ReS a (+i0) = because of the particle- hole symmetry, it is easy to see that 
1/R a (e + iO) has a pole at e = and that, therefore, S(e + iO) also has a pole at 
e = 0. HR a (e+iO) is given by Eq. (03]), it follows from Eq. gdj that ImS(£-N'0) = 
for |e| < e g /2 but e = 0. Thus, Eq. (|4-4j) can be described as 

S(e + iO) = — ^ + 5S(e + iO), (4-6) 

where \ s = -l/[(d/de)R a (e + iO)] > and 

8S(e + i0) = / dx + dx\ ^—, (4-7) 

V ' vr [J^ J +eg/2 J e + iO-x' V > 

where As(x) = A${— x) > 0. In order for a gap as large as e g to open in the Anderson 
model, S(e + iO) has to have a pole at e = and the residue Xs has to be so large 
that |e — Re5cr(e + i0)| cannot be zero for any e of |e| < e g /2; at leastHI 



Xs>~ 4 



el + 2e g \5S(e g /2 + iO)\], (4-8) 



has to be satisfied. 

Since S(e + iO) is the sum of S ff {e + iO) — /i and 7^(e + iO), as shown in Eq. (|4-2p . 
it follows that 

Z' (T (e + i0)-/i = ^- + ^ I :(e + i0), f(e + iO) = ^-+5S r (e + iO), (4-9) 
e + zL) e + iO 



where Xs + Ap = A5, and that 



1 / /' _eg//2 /" +00 

5Sz(e + i0) = — / dx+ dx 

1 / f~ e z/ 2 f +OQ 

5S r {e + i0) = — / dx + / dx 

71 \J-00 J+e s /2 



e + iO — x' 


(4-10) 


e + iO — x' 


(4-11) 



where Az(x) + ^lr(^) = As(x) for |x| > e^/2; in general, Xjj > 0, Xr > 0, /li:(x) = 
A~e{-x) > for |x| > e g /2, and A r {x) = A r (-x) = A(x) > for \x\ > e g /2. 

An issue is whether or not r(e + i0) can have a pole at e = 0, i.e., whether or not 
Ar can be nonzero. There are two possible cases for Xr'- Xr = and Xr > 0. First, 
we examine the case in which Xr = 0, i.e., X^j > and the single-site S ff (e + iO) has 
a pole at e = 0. According to the analysis in AppendixlA) Ar = means that the 



*' If the multisite AS a (e + iO, k) is continuous as a function of fc, Eq. (|3-8p or a similar one has 
to be satisfied for a complete gap to open in the Hubbard model, so that no complete gap can open 
if W is infinite. If it can be discontinuous as a function of k, a complete gap can open even if W is 
infinite, as examined for a continuous model in Appendix[C] If a complete gap opens in a finite or 
infinite bandwidth model, Eq. (|4-8[1 has to be satisfied in either model. 
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multisite AU a (e + iO, k) has no pole at e = 0. The possibility of this type of solution 
cannot be excluded within the qualitative analysis in this section. 

Second, we examine the case in which Xp > 0, i.e., A(e + iff) includes the delta 
function 5(e). In the original definition of r(e + iO), which is given by Eq. (|2-13a|) . 
we first assume that 



A(e) = \ r 



lr 



£ 2 + 7f 



(4-12) 



where jr is nonzero and finite such that < jp *C e g ; then we eventually take the 
limit of 7r —> 0. It is evident that Tk is neither zero nor infinitesimal if jr is neither 
zero nor infinite. If Tk defined by Eq. (|2-25|) is used, it follows from Eq. (|2-26|) that 



E <T (e + iff) = -U + 



l-nA(0)/(2W s k B T K ) 



e + 



(4-13) 



where W s is the Wilson ratio defined by Eq. (|2-30p . This expansion is relevant for 
\e\ < min(7f, k B Ta). It follows that 



lim 

7 r — 5-0 



de 



Re2J a (e + iff) 



lim 

7r— >0 



1 



7rAr 



2W s k B T Klr 



-co. 



(4- 14a) 



The Kondo temperature Tk is a function of jr- There are two possible cases for Tk 
in the limit of jr — > 0: /cbTk > 7r and /crTk < Jr- If ^bTk > Jr, then 



lim Re[i7 CT (-7r + iO) - K(j r + i0)] 

If /crTk < 7r, then 

lim Re\lJ a (-k B T K + »0) - £a(k B T K + »0)1 



TT~X r /(W s k B T K ) > 0. 



lim TrXr/iWsjr) 



(4- 14b) 



+oo. (4- 14c) 



According to Eq. (|4-14|> . the real part of [^(e + iO)] ~ Q is discontinuous at e = 0, 
so that its imagery part is divergent at e = in any possible solution of this type. 
If it continuously diverges as e — > such that 



lim|rmi7 CT (£ + i0) 



7r— »0 



-oo, 



(4-15) 



a zerogap opens, which contradicts the opening of a complete gap. A complete gap 
opens only if it discontinuously diverges at e = such that 



[imS^e + iO)] ac-S(e) 



(4-16) 



i.e., only if S a (e + iff) has a pole at e = 0. This is only possible if Tk — > OK in the 
limit of 7r — > 0. Since A(0) — > +oo as jr — > 0, i.e., A(0) > 0, it is probable that 
Tk = +0K rather than Tk = OK; at least, the entropy at T = +0K is nonzero. 
This type of solution is only possible beyond S 3 A, as studied in AppendixlAl 

In this section, we list all the possible types of solution for a complete-gapped 
phase beyond S 3 A in an unspecified model. It is shown that S ff (e + iff) has to have 
a pole at e = in any of them, which means that Tk = OK or Tk = +0K in it. 
Thus, the residual entropy or entropy at T = +0 K is nonzero in it; the entropy at 
T = +0 K is k B In 2 per unit cell, as studied in AppendixlBl 
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§5. Discussion 

Since the density of states p{e) is independent of wave number, it is essentially 
a single-site or local property in any model. It is therefore reasonably anticipated 
that p(e) can be described by that of a local or impurity model even for any model, 
even beyond S 3 A or DMFT; S 3 A and DMFT are exactly equivalent to each other. 
According to the Kondo-lattice theory™ >E^IIIJ,I23 the p(e) of the Hubbard model 
can be actually described by that of the impurity Anderson model, either under 
or beyond S 3 A, as shown in Eq. (|2-19bj) . Since the site-diagonal R a (iei) includes 
the total S a (iei,k), as shown in Eq. (J2-18J) . the mapping condition of Eq. (|2-14d|) 
depends on the multisite AE a (iei,k), i.e., the hybridization energy A(e) of the 
Anderson model to be mapped depends on AU a (isi, k); the multisite AU a (iei,k) 
is approximately or rigorously considered beyond S 3 A, although it is ignored under 
S 3 A. Beyond S 3 A, therefore, all of the A(s), £ a (iei), and A£ a (isi,k) have to be 
self-consistently calculated with each other to satisfy the mapping condition of Eq. 
(|2-14d|) . The effects of intersite correlations, which are described in terms of the 
multisite AU a (iei, k) of the Hubbard model, can be described in terms of the self- 
consistent A(e) of the Anderson model, at least, as regards any local property such 
as p(s). 

The residual entropy or entropy at T = +0 K is nonzero for any possible 
complete-gapped phase either under or beyond S 3 A or DMFT; in particular, the 
entropy at T = +0 K is k-g In 2 per unit cell, as studied in AppendixlBl In Brinkman 
and Rice's theoryp^* the metal-insulator transition appears continuous. According 
to the analysis in the present paper, it is discontinuous such that the residual entropy 
or entropy at T = +0K jumps at the transition. 

Theories based on DMFT^'^'ESli are consistent with Brinkman and Rice's 
theory except for the hysteresis in DMFT. Since the transition is discontinuous, it is 
reasonable that the hysteresis appears. Under DMFT, a metal-insulator transition 
occurs at U C2 as U increases, while it occurs at U Cl , which is smaller than U C2 , as 
U decreases. According to the analysis in the present paper, the hysteresis can be 
explained in the following way: the Kondo temperature Tk continuously vanishes at 
U C2 as U increases, while it discontinuously becomes nonzero at U Cl as U decreases; 
and the residual entropy or entropy at T = +0 K discontinuously becomes nonzero 
at U C2 as U increases, while it discontinuously becomes zero at U Cl as U decreases. 
The jump of Ik at U Cl means that the internal energy also jumps at U C1 . 

In general, an infinitely degenerate ground state is unstable even in the presence 
of an infinitesimal perturbation, except for a symmetry-broken state in which a 
continuous symmetry is broken and rigidity appears.^* This is why the third law 
of thermodynamics is valid. It is therefore doubtful that a gapped phase with no 
symmetry broken, in which the residual entropy or entropy at T = +0 K is nonzero, 
is stable under or beyond S 3 A in the presence of an infinitesimal perturbation due 
to the communication of electrons between the Hubbard model and an electron 
reservoir, i.e., in the grand canonical ensemble, as reported in a previous paper P» 

In two dimensions and higher, symmetry can be broken when the temperature 
T is low enough and, at least, at the absolute zero Kelvin. An infinitely degenerate 
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ground state in which the third law of thermodynamics is broken is unstable against 
an ordered state such as an antiferromagnetic state. The third law of thermodynam- 
ics is valid in any ordered ground state because rigidity appears in itP^J 

The Kondo temperature 7k or £;bTk is the energy scale of single-site quantum 
spin fluctuations. Electrons behave like local moments in a high-T phase at T ^> Tk, 
because single-site thermal spin fluctuations overcome single-site quantum ones. The 
high-T phase at T > Tk is in a sense a paramagnetic type of the Mott insulator, 
which is not the ground state. If an antiferromagnetic gap opens, the ground state 
can also be an insulator.^ A local-moment type of antiferromagnetic insulator, in 
which the Neel temperature Tn is as high as Tn S> Tk, is an antiferromagnetic type of 
the Mott insulator, because electrons behave like local moments in its paramagnetic 
phase at T > Tm . If electrons are not half- filled and the nesting of the Fermi surface 
is not sharp, Tn is low or it vanishes. If Tm is nonzero and Tn <C Tk, electrons 
are itinerant in a paramagnetic phase at Tn < T <C Tk, so that the magnetism at 
T < Tn is itinerant-electron magnetism. If Tk is high enough, it is possible that the 
ground state is superconducting^ 

The Hubbard model in one dimension is particular, because no symmetry can 
be broken in itp=* and because the Bethe-ansatz solution was given by Lieb and 
Wu.Hs3 \Ye denote the ground-state energy in the canonical ensemble by Eq(N) as 
a function of the number N of electrons. Two chemical potentials are defined for the 
ground state of N electrons in the canonical ensemble: /i+(iV) = Eq(N+1) — Eq(N) 
for the addition of an electron and /U_(iV) = Eq(N) — Eq(N — 1) for the removal 
of an electron. A gap is defined by e s (N) = fi + (N) — fi_(N). According to the 
Bethe-ansatz solution, e g (iV) = for N ^ L, while e g (-/V) > for N = L for any 
nonzero U; i.e., a gap opens in the spectrum of adding or removing the whole of 
a single electron in the exactly half-filled case. On the other hand, the third law 
of thermodynamics is not broken in the ground state of the Bethe-ansatz solution. 
If the opening of a gap in the spectrum of adding or removing an electron in the 
canonical ensemble is equivalent to that in the single-particle excitation spectrum in 
the grand canonical ensemble, the conclusion of the present paper and the argument 
based on the Bethe-ansatz solution are contradictory to each other. 

A possible explanation for the contradiction is that the Kondo-lattice theory, 
which is none other than the 1/D expansion theory from infinite dimensions, is 
inapplicable to the opposite limit of one dimension. Since the Kondo-lattice theory 
is rigorous in infinite dimensions as S 3 A is, it is applicable at least to D — > +oo 
dimensions^ If the explanation is right to the point, therefore, there is a critical 
D c such that the Kondo-lattice theory is applicable to D > D c dimensions while it 
is inapplicable to 1 < D < D c dimensions. It is desirable to determine how large or 
small the critical D c is. It is also desirable to clarify why the Kondo-lattice theory is 
inapplicable to 1 < D < D c dimensions and, in particular, why the mapping of the 
single-site self-energy of the Hubbard model to the local self-energy of the Anderson 



*' If the multisite A£ a {e+iO, k) as well as the single-site S a (e + iO) is rigorously calculated, the 
Kondo-lattice theory that includes no conventional Weiss mean field, if it is applicable and relevant, 
is rigorous for any D dimensions, but within the constraint Hilbert subspace in which no symmetry 
is allowed to be broken the same as S 3 A for infinite dimensions. 
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model is impossible for 1 < D < D c . Even if the Kondo-lattice theory is really 
inapplicable to 1 < D < D c dimensions, a necessary condition for the opening of a 
gap in the single-particle excitation spectrum is that the total self-energy U a (e+iO, k) 
has a pole at e = at least for k = k-p on the hypersurface defined by Eq. (|2-4p . as 
discussed in the last paragraph in £ )2.2i Since the self-energy E a {e + iO,k) having 
a pole at e = implies or at least suggests that there are degenerate zero-energy 
bosonic excitations, which are responsible for the divergent imaginary part of the 
self-energy, it is desirable to elucidate why and how the residual entropy can be zero 
even in such an anomalous situation. 

Another possible explanation for the contradiction is that the opening of a gap 
in the canonical ensemble is not equivalent to that in the single-particle excitation 
spectrum in the grand canonical ensemble. Although the above well-known argu- 
ment by Lieb and Wu is certainly based on the exact result of [eg(N)] N=L > 0, 
which is obtained by the Bethe-ansatz solution, the inequality of e g (N) > is not 
a sufficient condition for the ground state being an insulator; for example, the very 
same inequality is satisfied in a metallic fine particle if the long-range Coulomb inter- 
action is explicitly considered. Thus, the possibility that the half-filled ground state 
is a metal is not excluded by the Bethe-ansatz solution itself. It should be examined 
whether a gap opens in the single-particle excitation spectrum in the grand canonical 
ensemble. If the Kondo-lattice theory is applicable to one dimension, a combination 
of the ground-state degeneracy of the Bethe-ansatz solution and the conclusion of 
the present paper means that no complete gap opens in the single-particle excitation 
spectrum. It is desirable to study the conductivity itself in the presence of electrodes, 
which plays a role of an electron reservoir for the Hubbard model. 

In the continuous limit of a — > 0, where a is the lattice constant, the Hubbard 
model in one dimension is substituted for the Tomonaga-Luttinger model with umk- 
lapp terms^'^'^ and the Tomonaga-Luttinger model is further substituted for a 
boson model, in which a charge gap opens. 461 Since the Tomonaga-Luttinger model 
is a continuous or nonlattice model, the Kondo-lattice theory is inapplicable to it, 
i.e., the Kondo effect is irrelevant in it, as discussed in AppendixICl If the Kondo- 
lattice theory is applicable to the Hubbard model in one dimension, the opening of a 
complete gap in the Tomonaga-Luttinger or Boson model, in which the Kondo effect 
is irrelevant, is not a proof that contradicts the conclusion of the present paper for 
the Hubbard model, in which the Kondo effect is relevant. 

§6. Conclusion 

The self-energy of the Hubbard model is decomposed into the single-site and 
multisite ones. According to the Kondo-lattice theory, the single-site self-energy is 
equal to that of the mapped Anderson model that is self-consistently determined 
in a way such that the density of states of the Hubbard model per unit cell is 
equal to that of the Anderson model. If the multisite self-energy is ignored under 
the supreme single-site approximation (S 3 A) or if it is approximately or rigorously 
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considered beyond S 3 A, the density of states per unit cell is simply given by 

p(e) = Im = , (6-1) 

where S(s+i0) = e^ — p,+E a (£+i0)+r(£+i0), where id, jl, £ a (£+iQ), and T(e + i0) 
are the level of localized electrons, the chemical potential, the self-energy, and the 
hybridization term, respectively, of the mapped Anderson model; and fi—id = f J -~ e d, 
where \i and e^ are the chemical potential and band center, respectively, of the 
Hubbard model. The Kondo temperature Tk or &bTk of the mapped Anderson 
model is also the low-energy scale of single-site quantum spin fluctuations in the 
Hubbard model. 

In the exactly half-filled case of the Hubbard model on a bipartite lattice, 
Re5 r (+i0) = because of the particle-hole symmetry. A complete gap as large 
as e g can open in p(e), if and only if S(e + iO) has a pole at e = 0, the residue of 
the pole is so large that \e — ReS a (e + iQ)\ cannot be zero for any |e| < e g /2, and 
ImS(s + iO) = for any |e| < e g /2 but e = 0. Such an anomalous S(s + i0) is possible 
if and only if Tk = K or Tk = +0 K. If Tk = K, the residual entropy is nonzero. If 
Tk = +0K, the ground state is a singlet, while the entropy at T = +0K is nonzero; 
the entropy discontinuously vanishes at T = K with decreasing T. Thus, the Mott 
insulator with no symmetry broken, if it is possible, is characterized by nonzero 
residual entropy or nonzero entropy at T = +0 K. 

The conclusion of the present paper is consistent with Brinkman and Rice's 
theory, in which the specific-heat coefficient diverges as the onsite U approaches a 
transition point from a metallic phase and no symmetry is broken in an insulating 
phase; a combination of these two properties means that the residual entropy or 
entropy at T = +0 K is nonzero in the insulating phase. It is also consistent with 
the appearance of hysteresis in the dynamical mean-field theory. 

According to the well-known argument based on the Bethe-ansatz solution, on 
the other hand, the half-filled ground state in one dimension is the Mott insulator 
for any nonzero onsite U, although the third law of thermodynamics is not broken in 
it; thus, the well-known argument and the conclusion of the present paper contradict 
each other. A possible explanation for the contradiction is that the Kondo-lattice 
theory, which is none other than the 1/D expansion theory from infinite dimensions, 
is inapplicable to one dimension. The other possible explanation is that a gap in 
the spectrum of adding or removing an electron in the canonical ensemble, which 
is given by the Bethe-ansatz solution, is different from a gap in the single-particle 
excitation spectrum in the grand canonical ensemble, which is studied in the present 
paper. Since the opening of the gap in the canonical ensemble is not a sufficient 
condition for the ground state being an insulator, the well-known argument is not a 
proof that the half-filled ground state in one dimension is the Mott insulator. If the 
Kondo-lattice theory is applicable to one dimension, a combination of the conclusion 
of the present paper and the ground-state degeneracy of the Bethe-ansatz solution 
means that the half-filled ground state in one dimension is not the Mott insulator. 
It is desirable to elucidate which explanation is true in the near future. 
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Appendix A 

Anomalies of Gapped Phases 



In order for a complete gap as large as e g to open, the single-site E a {e + iO) has 
to have a pole at e = 0, as examined in §§ 13.21 and [U so that the total U a {e + iO, k) 
also has to have a pole at e = 0. However, it does not have to have two or more 
poles in the gap region of |e| < e g /2, as discussed below. 

If the total 2J a (e + iO, k) has two or more poles in the region of \e\ < e g /2, there 
is at least a pair of adjacent poles: e.g., one at e = p\ and the other at e = p2, where 
— e g /2 < pi < P2 < £g/2- Then, it follows that 

lim ReS a (e + iO, k) = +oo, lim ReU a (e + i0,k) = -oo. (A-l) 

£— >Pl+0 £— >P2 — 

Since U a (e + i0, k) has no singularity between poles and is continuous between poles, 
there is a zero point z £ between poles such that it satisfies 

Re[z £ + ijl - E{k) - U a (z £ + i0,k)] =0. (A-2) 

Since the zero point z e is in the gap region of |e| < e g /2, it follows that 

lm[z £ + fi- E{k) - U a (z £ + iO, k)] = 0. (A-3) 

Since the density of states p(e) of the Hubbard model is given by Eq. (|2-19p with 
Eq. (|2-18p . an in-gap state inevitably appears at each zero point e = z £ between each 
pair of adjacent poles. This conclusion contradicts the opening of a complete gap. 
Thus, the total self-energy E ff {e + iO, k) does not have to have two or more poles in 
the gap region of |e| < e g /2 for a complete gap as large as e g to open. 

A complete gap as large as e g opens in the Hubbard model if and only if single- 
site and multisite self-energies are given by 



2J a (iei) - n= -. 

lEl IT 



AS a (iei,k) 



\Az(k) 



i£l 




(A-4a) 

(A-4b) 



where Xs + ^Ai:(k) is nonzero and large enough; \% > and \As(k) > 0, in general. 
According to Eqs. (pMKD and (lA4l> . 



RJieA 



lEl 



\Y. 



+ 



faY 



(A-5) 



in the limit of |e/| — > 0. According to the mapping condition of Eq. ()2-14d|) . 

A(e) = (XrM5(e), (A-6) 
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for |g| < e g /2, where 

IE 



X r 



1 



1 



Aj; + AAc(fc) A^ 






Ar > 0, A^(fc) > 0, 



\e = 0, X A z{k) > 0. 



(A-7) 



If \Az(k) = or AU a (e + iO,k) has no pole at e = 0, then Xr = 0, i.e., /"(e + iO) has 
no pole at e = 0. If A/ix^fc) > or AS a (e + zO, fe) has a pole at e = 0, then Xr > 0, 
i.e., r(e + iO) has a pole at e = 0. According to the analysis in £JU if r(e + iO) 
has a pole at e = 0, the single-site E a {e + iO) also has a pole at e = 0. Thus, any 
self-consistent solution characterized by Xs = and XAs(k) > is impossible. 

Next, we assume that the single-site E a {e + zO) is divergent at e = such that 

]im.\£Je -HO) I = +oo, \SJ+i0)\ = +oo, (A-8a) 

while the multisite AE a [e + iO, k) is finite and continuous at e = and 

Im[Z\E CT (e-H0,fc)] =0, (A-8b) 

for any |e| < e g /2 and any fc. We define 



9(e) 



e + i0 + fi- £ a (e + i0) 
It is easy to see that 



Qn{e) = jYl i E W + AZ a {e + iO, k)] n . (A-9) 



lim|$(e)|=0, |g(0)|=0, 



(A-10) 



Q n (e) is real for |e| < e g /2, and Q n (e 
Eq.CZHBJ, 



(— l) n Q n (— e) for |e| < e g /2. According to 

fitj(g X + . Q) = ~^y " Qi(e) " S(e) [Q2OO - Q?00] + O [^(e)] , (A-ll) 

for a sufficiently small |e|. According to the mapping condition of Eq. (|2T4d|h 

A(e) = -Im{g(e)[Q 2 (e) - Qj(e)] + 0[g\e)}}, (A-12) 

for a sufficiently small |e|. It follows from Eqs. ()A-10p and (|AT2p that 

limA(e) = 0, A(0) = 0. (A- 13) 

In particular, the multisite AE a (e + iQ, k) vanishes under S 3 A. Therefore, Eq. (|AT3p 
has to be satisfied in any possible complete-gapped or zerogapped phase under S 3 A. 
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Appendix B 

Entropy at T = +0 K of Possible Complete- Gapped Phases 



If no symmetry is broken, the thermodynamic potential is given by-^'E^ 1 
/? = -fc B r^^e^ 0+ {ln[-l/G (7 (^,/c)] +G <T (i£ h k)Z <7 (ie l ,k)} + f2', (B-l) 



l kcr 



where Q 1 is the contribution of all the closed connected Feynman diagrams but with 
each electron line replaced by its renormalized one and is given by 

Q' = k B T Y^ E E ^ G ^' k ) E -^ i£ h k), (B-2) 

v I kcr 

where E ua {iei,k) is the z/th order self-energy in U, where only the U occurring in 
explicit interactions of the skeleton diagram is used to determine the order. Since 
the thermodynamic potential Q is stationary with respect to the self-energy such 
that df2/dU a (iei,k) = in the functional derivative, the temperature dependence 
of U a (iei,k) does not have to be considered in calculating entropy S = —dQ/dT 
except for that coming from iei or the I sum. 

As discussed by Luttinger and Wardp^'™ the first correction to the T = value 
of Q comes only from the difference between the I sum and what we would obtain if 
we replaced them by integrals in each electron line. Then, the first correction to the 
T = value of fl' is equal to that of 



n> = k B TY,Y,Y, G ^ i£ ^ k ^^ i£ ^ k ^ 

v I kcr 

= k B TY,Y. G °( i£l >^ I] Ai£i,k), (B-3) 



l kcr 

so that the first correction to the T = value of the total Q is equal to that of 
Q{T) = -feTEE^W-VGU^fc)] 



l kcr 

l-OO 



i r+oo 

= ~E/ def{e)lmln[-e-iO + E(k) + U a {e + iO,k)-fi], (B-4) 

^ ka J -°° 

where /(e) = l/[ e 6 /^ T ^ + l] ; in Eqs. (|B~3l) and (TB4l) . T = OK is assumed in the 
self-energy except for the / sum. 

For any possible complete-gap phase either under or beyond S 3 A, it follows 
according to Eq. (|A-4j) that 

SAe + „, k) _ , = MS _ i ( n l2 dx + r dx \ m^l, (B . 5) 

e + iO vr yj_ 00 J +eg/2 J E + lO-X 

where \s(k) = Xe + A\ A z(k) > and A E (x; k) = A s + AA A s(x; k) > 0. Then, it 



o, 


-Cg/2 < e < 0, 


vr/2, 


e = 0, 


-7T, 


< e < +e g /2, 
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follows that 

Imln[-e-*0 + J5(fe) + i7 ff (eH-«0,fe)-/x] = < -vr/2, : = 0, (B-6) 

for |e| < Cg/2. The entropy at T = +0K per unit cell is as large as 

S/L = - lim -^-1= lim ^- / def{e) = k B \n2, (B-7) 

for any possible complete-gap phase either under or beyond S 3 A. 

Appendix C 

Continuous Model in One Dimension 



We consider a finite or infinite bandwidth model in one dimension in the con- 
tinuous limit of a — > 0, with any of k?a = tt/2, hv-p = \{d/dk)E{k)\ fc=fc , g = Ua, 
and La being kept constant; vp is the Fermi velocity and g is the coupling constant. 
We assume the thermodynamic limit of La — > +oo. 

If a is nonzero and finite, the single-site Green function is given by 



a r+ n / a 1 

RJe + iO) = — dk — — — , (C-l) 

V ; ^J-Kla e + iO-E(k)-S(e + iO,k)' v ' 



(tJ r'KJa -1 

'— 7r/a 

where 

S(e + i0,k) = U a (e + iO,k) - //. (C-2) 

Since it vanishes in the continuous limit of a — > 0, the single-site self-energy is never 
well-defined in any continuous model. Thus, the Kondo-lattice theory is inapplicable 
to any continuous model, so that the Kondo effect is irrelevant in it. 

When the infinite bandwidth model defined in ^2.11 is considered, for example, 
the bare density of states per unit length is given bvT-1 

Po(E) = - P0 (E) = - F2 ^°' = I t-tttAtz , 12 , (C-3) 

a ira E z + \t QO \ jL Trnvp 1 + [Ea/(hvY)\ 

where a = +0. Since E(k) is a continuous function of k, k is also a continuous 
function of E. Thus, S[e + i0,k(E)] is a function of E; it is simply denoted by 
S(e + 10, E). The density of state per unit length is given by 

p(e) = -J^ ^(E) E + i0 _ E _ s{s+ . 0Ey (04) 

Because of the particle-hole symmetry, ReS(-M0, 0) = 0. It is easy to see that unless 
ImS(e + iO, 0) is divergent at s = 0, then p(0) > and not p(0) = 0. 



* 5 If po(E) = e~ alEl/{hVF) /(Trhv F ) with a = +0 is assumed instead of Eq. (fCT3)) , it is none other 
than the band cutoff introduced into the boson model. 



Residual Entropy of the Mott Insulator with No Symmetry Broken 



25 



We examine a possible complete- gapped phase whose gap is as large as e g : p(e) = 
for |e| < e g /2. S(e + iO, E) as a function of E may be continuous or discontinuous. 
First, we assume that S(e + iO, E) is continuous as a function of E. According to 
an analysis similar to that in AppendixlAl S(e + iO, E) does not have to have two or 
more poles in the region of |e| < e g /2 for a complete gap to open. Thus, it has to 
have zero or only one pole in the region, so that 



S(e + iO,E) 



X(E) 



A(x;E) 



e + iO-p(E) vr 




e + iO 



(C-5) 



where X(E) = X(-E) > 0, p(E) = -p(-E), and A(x;E) = A(-x;-E) > 0; in 
particular, A(0) > and p(0) = 0. If the pole p(E) depends on E, S{e + iO, E) as a 
function of E is discontinuous at E such that it satisfies p(E) = e. Thus, p(E) = 
because of the constraint imposed by the assumption above. It is evident that 



lim [E + S(e + iO,E)] = +oo, lim [E + S(e + iO, E)] 

E— >>+oo ' E— >— oo 

Then, there is a zero point ze such that it satisfies 

Re [e- z E - S(e + iO, z E )} = 0, 



-oo. 



(C-6) 



(C-7) 



for any finite e but e = 0. The zero point ze is a function of e; the function is 
denoted by ze(s). If |e| < e g /2 but e 7^ 0, 



Im{e- z E {e) - S [e + i0,z E (e)]} = 0. 



(C-8) 



Then, p(e) > for any e but e = 0. This conclusion contradicts the opening of 
a complete gap. If the bandwidth is infinite, no complete gap can open in p(e); 
however, a possibility that a zerogap opens in p(e) cannot be excluded. 

If S(e + iO, E) is allowed to be discontinuous as a function of E, numerous types 
of anomaly are possible because S(e + iO, E) as a function of E does not need to be 
analytical in the upper or lower half complex plane of E; for example, 



S(e + i0,E) 



X(E) 



e + iO 
where X(E) > and, in particular, A(0) > and 



+ 5S(e + i0,E), 



(C-9a) 



5S(e + iO, E) 



E 

Te\ 



E 2 + (eg/2) 2 _ |^| 
0, 



E^0, 

E = 0, 



(C-9b) 



S(e + i0,E) 



(%/2) 2 



e + iO + E ' 



(C-10) 



and so on. If S(e + iO, E) is given by Eq. (|C-9p or (|C-10p . a complete gap as large as 
e g opens, even if the bandwidth is infinite. A possibility that a complete gap opens 
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in p(e) cannot be excluded. A necessary condition for a complete gap to open is that 
S(e + iO, E = 0) or S(s + i0,k = ±kp) has a pole at e = and S [e + iO, E(k)] is 
discontinuous as a function of i? or k. 

The single-band model examined above can be substituted for a two-band model. 
We define two types of creation and annihilation operators by 



a k 



ka ~ d (k+k F )a> a kT ~ d(k+k F )a, P k a - ^(fc-fep)^' Pkff - d(k-k F )a, (C-ll) 



for — kp < k < ky; a and (3 electrons correspond to right-going and left-going ones, 
respectively, in the single-band model. The one-body term is substituted for 



H = Y, M («L«far " Pljke) ■ (C-12) 

ka 



The two-body term is decomposed into several terms according to the numbers 
of a k(T , aka, b ka , and bk a , some of which are the so called umklapp terms. This 
two-band model is none other than the Tomonaga-Luttinger model with umklapp 
terms.'^Z'"^'''^ Since the two-band model is a continuous model the same as the 
original single-band model, the Kondo-lattice theory is inapplicable to it. Since each 
of the numbers of a and /3 electrons is not a conserved quantity because of the 
two-body term, the Green function is given, in general, by 

q i ie k \ = ( isi-fvopk- S aaa {ie h k) -S al3a {ie h k) \ (013) 

''"'" V Sfiaff(i£i,k) isi + hv F k- Si3i3 a (iEi,k) J 

If the correspondence to the original single-band model is seriously considered, the 
interband terms have to vanish such that 

S a p ff (iei,k) = Sp aa (i£ h k) = 0, (014a) 

and the intraband terms have to satisfy 

S aaff {i£i,k) = Spp ff (i£i,-k) = £ a [i£i,E(k + k F )] - fi. (014b) 

The density of states of this two-band model is equal to that of the original single- 
band model, provided that Eq. (jQ14j) is satisfied. Thus, a necessary condition for a 
complete gap to open is that S aaa {e + iO, k) has a pole at e = at least for k = 
and S aaa (e + iO, k) is discontinuous as a function of k. 

It is desirable to elucidate how discontinuous S aaa (s + iO, k) is as a function of 
k in the gapped phase of the Tomonaga-Luttinger model, which corresponds to the 
gapped phase of the boson model. 46 ' Since it is surprising that the self-energy as 
a function of wave number k is discontinuous, it is desirable to elucidate what is 
actually responsible for the discontinuity. It is also desirable to elucidate whether 
the residual entropy is zero or nonzero in the gapped phases; the divergence of 
ImiS r aao -(H-iO, 0) implies or at least suggests that there are degenerate zero-energy 
bosonic excitations, which are responsible for the divergent ImS aaa (+iO,0). 
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